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The  combination  of  time  domain  source  imaging  techniques  and  arrays  to  consider  the  localization 
of  an  impulsive  source  in  a  wedge  waveguide,  with  a  free  surface  and  rigid  bottom,  is  investigated. 

An  impulsive  signal  is  transmitted  from  an  unknown  location,  received  by  an  array  of  hydrophones, 
and  stored.  The  signals  are  time  reversed  and  transmitted  into  the  model  wedge  by  a  transducer  at 
each  receiver  location.  Model  transmission  functions  are  used  to  compute  the  pressures  at  a  set  of 
locations.  The  field  calculations  use  an  image  construction  from  D.  Chu’s  exact  solution  for  a 
density  contrast  wedge  [J.  Acoust.  Soc.  Am.  86,  1883-1896  (1989)].  A  map  of  the  peak  pressures 
gives  an  image  of  the  source  location.  Particular  attention  is  given  to  the  peak  amplitudes,  sidelobe 
amplitudes,  and  spatial  resolution  as  a  function  of  the  number  and  placement  of  receivers.  Source 
localization  is  enhanced  by  the  range  dependency  of  the  environment,  which  eliminates  “range” 
sidelobes,  even  for  a  single  hydrophone.  Adding  receivers  along  an  arc  eliminates  “angle” 
sidelobes;  but  adding  receivers  along  an  arc  or  radial  gives  little  reduction  in  overall  background 
level.  Adding  receivers  parallel  to  the  wedge  axis  gives  sharp,  unambiguous  azimuthal  imaging. 

PACS  numbers:  43.20.Mv,  43.30. Vh,  43.30. Wi,  43.35.Sx 


INTRODUCTION 

There  are  two  extreme  types  of  sources:  Sources  of  con¬ 
tinuous  wave  harmonic  pressures  and  sources  of  broad  fre¬ 
quency  bandwidth  pressures.  Source  location  methods  differ 
for  the  two  kinds  of  sources.  Matched  field  techniques  for 
harmonic  sources  require  an  array  of  receivers.1-4  Matched 
signal  methods  use  wide  frequency  bandwidth  signals  and  as 
few  as  a  single  receiver.5-10  Both  methods  use  many  obser¬ 
vations  of  the  field  in  the  waveguide  but  do  it  differently. 
Matched  field  methods  typically  use  the  observations  at  an 
array  of  receivers  and  a  single  frequency.  Matched  signal 
methods  may  use  observations  at  many  frequencies  and  one 
receiver.  Both  methods  give  ambiguity  or  image  maps  and 
these  maps  can  have  large  “false”  images  or  peaks  that  rival 
or  dominate  the  proper  image  at  the  source  location. 

Since  matched  signal  methods  can  use  a  single  receiver 
to  get  a  source  location,  does  an  array  of  receivers  “im¬ 
prove”  the  image  or  ambiguity  maps?  There  are  many  ways 
to  measure  any  improvement  that  may  be  obtained.  Here  we 
choose  to  use  the  ratio  of  the  maximum  of  the  false  (or 
“sidelobe”)  peaks  to  the  peak  at  the  matched  position  or 
source  location.  The  spatial  resolution  of  the  image  is  also  of 
interest.  Numerical  studies  of  source  location  in  a  Pekeris 
waveguide  have  shown  that  the  false  peaks  were  reduced  by 
increasing  the  number  of  receivers  in  an  array;  however,  the 
results  were  not  dramatic.10  Matched-mode  methods  and 
time  domain  methods  can  be  combined.11 

As  an  extension  of  our  previous  research  in  a  waveguide 
having  parallel  boundaries,10  we  study  the  wedge  waveguide. 
To  emphasize  the  role  of  the  receiver  array  in  the  image 
mapping  process  and  simplify  the  numerical  work,  we  use  an 
ideal  wedge  waveguide.12  The  exact  solutions  for  the  wedge 
do  not  require  intensive  numerical  computations.  For  the 


same  reason,  we  use  a  pulse  transmission  instead  of  a  ran¬ 
dom  source  with  the  same  signal  bandwidth. 

The  acoustical  source  location  method  used  here  is,  us¬ 
ing  an  optical  analogy,  a  holographic  imaging  technique.  The 
array  of  receivers  samples  and  records  the  outward  traveling 
wave  field  of  the  source.  If  recorded  signals  are  transmitted 
by  an  array  of  transducers  at  the  receiver  locations,  then  the 
wave  field  continues  traveling  outward  as  if  it  came  from  the 
source.  If  the  recorded  signals  are  time  reversed  and  then 
transmitted  by  the  array  of  transducers,  the  backward  travel¬ 
ing  wave  fields  converge  to  form  an  image  of  the  source  at 
the  source  location.  The  time  reversal  technique  is  also 
known  as  phase  conjugation  because  the  complex  conjugate 
of  the  spectral  function  gives  a  time  reversal  of  the 
signal.12'14  Imaging  techniques  are  highly  developed  in  ex¬ 
ploration  seismology  where  they  are  known  as  “migration”, 
“imaging,”  and  “impulse  holography”.1315  Ideally,  the  back 
projected  wave  fields  are  shown  as  a  moving  picture  as  the 
waves  travel  in  the  waveguide.  Instead  of  showing  a  motion 
picture  or  a  sequence  of  snapshots,  we  make  a  quasistatic 
image  of  the  image  field  by  displaying  the  maximum  abso¬ 
lute  values  of  the  wave  fields  at  a  closely  spaced  grid  of 
points. 

Exact  solutions  for  wave  fields  in  a  wedge  are  in  Biot 
and  Tolstoy16  and  Chu.17  While  this  theory  provides  a  basis 
for  doing  the  required  calculations  exactly,  the  computations 
are  fast  and  easily  programmed.  The  method  facilitates  the 
investigation  of  source  imaging  in  the  wedge  as  a  function  of 
bandwidth  (by  varying  the  time  duration  of  the  transmitted 
pulse)  and  the  number  and  placement  of  receivers  (by  per¬ 
forming  reciprocal  calculations  to  multiple  receivers  at  dif¬ 
ferent  locations  in  the  wedge  and  aggregating  the  results). 
While  we  use  the  formulation  of  Biot  and  Tolstoy  and  Chu  in 
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this  paper  because  it  is  simple  and  direct,  Buckingham18  has 
given  a  different  formulation  for  calculating  the  sound  field 
in  a  wedge  due  to  continuous  wave  source^  !,;s  method 
included  reflection  beyond  the  critical  angle,  i  e  beam 

displacement  technique  described  by  Brekl  V''  and 
used  in  several  applications.20  2< 

Section  1  of  this  paper  gives  a  brief  review  of  useful 
expressions  from  wedge  theory,  image  construction  rays,  and 
the  time  reversed  transmission  technique.  Section  II  dis¬ 
cusses  the  imaging  of  acoustic  fields.  Section  III  describes 
the  environmental  model  and  signal  transmissions.  Section 
IV  describes  measures  of  localization.  Section  V  discusses 
single  receiver  and  pulse  duration  effects.  Section  VI  dis¬ 
cusses  array  designs.  This  is  followed  by  a  summary  of  the 
conclusions  drawn  from  the  study.  We  have  included  an  ap¬ 
pendix  which  describes  the  geometrical-kaleidoscopic 
method  for  constructing  reflected  ray  paths.  This  shows  how 
ray  path  calculations  may  be  performed  in  wedges  with  ar¬ 
bitrary  reflection  coefficients  at  the  boundaries  and  also  in 
finite  wedges. 

I.  THEORY:  ACOUSTIC  FIELDS  IN  A  WEDGE 

The  time  domain  impulsive  solution  for  the  wedge  prob¬ 
lem  has  tremendous  advantages  compared  to  the  single  fre¬ 
quency  harmonic  normal  mode  solution.  The  Biot-Tolstoy 
time  domain  expressions  are  a  sophisticated  solution  of  the 
wave  equation.1"161724  In  the  time  domain,  all  of  the  reflec¬ 
tion  arrivals  and  the  diffraction  from  the  wedge  apex  are 
identifiable  and  separated.  Depending  on  the  practical  situa¬ 
tion,  we  can  keep,  modify,  and  discard  any  of  the  arrivals. 

The  theory  for  the  wedge  was  given  as  an  example  of 
normal  coordinate  applications  to  problems  in  wave  propa¬ 
gation  [Biot  and  Tolstoy  (1957)].16  The  normal  coordinate 
method  comes  from  analytic  mechanics  and  is  often  called 
the  Lagrangian  method.  Readers  are  referred  to  Chap.  3  and 
Appendix  5  of  Tolstoy  and  Clay,  Ocean  Acoustics  (1987), 12 
and  Chap.  8  of  Tolstoy,  Wave  Propagation  (1973)24  for  nor¬ 
mal  coordinate  methods  in  wave  propagation.  The  original 
Biot  and  Tolstoy  solution  was  for  an  impulsive  point  source 
in  a  wedge  that  has  rigid  boundaries.  Later,  Kinney  e/a/.25 
gave  a  solution  for  the  wedge  with  free  boundaries.  Chu 
derived  the  exact  impulse  response  of  a  density  contrast 
wedge.17  He  also  derived  expressions  for  a  half-space  of  air 
over  a  half-space  that  contains  a  density  contrast  wedge.26 
This  paper  also  contains  an  extensive  list  of  references  on  the 
wedge  problem.  For  the  reader’s  convenience  we  give  here  a 
synopsis  of  solutions  for  the  different  boundary  conditions. 
In  the  shallow  water  waveguide  and  similar  problems  the 
diffraction  arrival  from  the  wedge  apex  is  likely  to  be  unim¬ 
portant  relative  to  the  contributions  from  the  image  part  of 
the  solution. 

The  coordinates  and  geometry  are  shown  in  Fig.  1.  The 
wedge  apex  is  along  the  y  axis.  The  boundary  planes  of  the 
wedge  are  normal  to  the  x-z  plane.  The  angles  are  measured 
counter  clockwise  from  the  lower  wedge  face.  These  choices 
let  us  use  the  z  coordinate  for  depth  and  to  have  the  x  and  y 
coordinates  in  the  horizontal  plane.  The  source  is  in  the  x-z 
plane  at  y=0.  The  receiver  can  be  displaced  a  distance  y. 
The  radial  distance  from  the  y  axis  is  r.  The  wedge  angle  ©*, 


FIG.  I.  Geometry  and  coordinates  for  wedge.  In  this  figure  the  coordinates 
of  the  source  are  labeled  and  those  of  the  receiver  (r.v.'Ft,  where 

the  angles  'F„,  'F  are  measured  down  from  the  top  (free)  surface  of  the 
wedge.  The  wedge  angle  is  B* .  The  angles  B„,  B  (for  the  source  and 
receiver,  respectively)  used  in  Sec.  1  are  measured  up  from  the  bottom 
(rigid)  surface  of  the  wedge.  They  are  B„=Bu  -'F()  and  B  =  BW  -'F. 

can  h —  any  value,  O<0m,«27t.  The  sound  velocity  is  c. 

■;t  and  Tolstoy1216  and  Tolstoy24  use  the  dis¬ 
play  ,.i  potential  </>  and  solve  for  r)<f>/r)t  for  a  unit  step  of 
displaces  ,  „t  the  origin.  Medwin  replaced  the  step  func¬ 
tion  in  the  volu'.ie  source  by  a  source  that  starts  a  constant 
flow  S  at  t=0.27  Then,  the  pressure  response  p(t)  is 

p(t)  =  {[P()  Ru]f>U-R/c)}/R  (1) 

for  pressure  P„  at  range  /?„.  We  write  hese  as  [P0  /?„] 
because  they  represent  the  source  pressure  and  always  go 
together  to  keep  the  dimensions  correct.  Th_  change  from 
d<t>/dt  to  pressure  from  an  impulsive  pressure  soiree  uses 

dtp 

/></)=— 4tt[P0  P()]  —  .  (2) 

The  derivation  of  the  complete  solution  of  the  wedge 
problem  is  long.  To  facilitate  the  reader’s  reference  to  the 
Biot-Tolstoy  and  Tolstoy  developments  we  use  d<f>/dt  and 
then  convert  to  p(t)  later.  The  changes  from  Biot-Tolstoy  to 
the  notations  in  Tolstoy,  Chu,  and  here,  are  a— >c,  /?— >y, 
z— >y,  f— >0* ,  and  x— >£.  Three  regions  of  the  solution  are 
defined  by  the  values  of  /  relative  to  t0  and  r0 

/0  =  (l/c)[(r-r0)2+y2]l/2,  (3) 

r0  =  (l/c)[(r  +  r())2  +  y2]1/2.  (4) 

We  define  the  fields  in  each  “time”  region  as 

<?<£  _d<t>  1  d<t>2 

dt  ~  at  +  dt  +  dt  ’ 

where  the  regions  occur  in  the  time  intervals  /</0, 
/„=£/<  r0,  and  T(j<f . 

A.  Region  1:  Before  the  first  arrival  and  t<t0 

Region  1  needs  little  discussion  because  <?</>,/<?/  is  zero 
prior  to  the  arrival  of  a  disturbance  from  the  source.  This  is 
equivalent  to  requiring  that  the  solution  be  causal 

d<f>  | 

~  =0,  /</„.  (6) 

B.  Region  2:  Direct  arrival  and  images  and  f„«f<T0 

Here  the  solution  is  the  sum  of  a  finite  number  of  im¬ 
ages.  The  expression  for  fields  in  a  wedge  having  rigid 
boundaries  follows: 
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FIG.  2.  Image  construction  for  a  free-rigid  shallow  water  wedge  waveguide. 
Wedge  angle=ll  deg.  Source:  angle  8  deg  from  surface  and  3000  m  range, 
y=0.  Receiver:  4  deg  from  surface  and  1000  m  range,  y=0. 


where  ct  becomes  Rk(+-)  etc.  and  these  are  given  by 
Rk(  +  -)  =  [rl  +  r2  +  y2-2r0r  cos(2*©h,+  0 - 0O)]I/2, 
/?*(-  +  )  =  [r^  +  r2+y2-2r0r  cos(2Ar©K,-0  +  0o)]1/2, 

Rk(  +  +)  =  [rl  +  r2+y2-2r0r  cos(2k@w+  0  +  0O)]I/2, 
Rk(~  ~)  =  [rl  +  r2  +  y2-2r0r  cos(2Ar©H,-@-©0)]1/2. 

The  maximum  value  of  k  or  A:max  is  given  by  the  condition 

O«|2it©H,±0±©o|:Srr.  (8) 

This  condition  gives  a  finite  number  of  images. 

If  one  lets  k=  -A:max,...0,...A:max,  then  (2A0W+0O— 0) 
and  (2^0*,— 0O— 0)  give  all  of  the  combinations.  The  pres¬ 
sure  disturbance  from  the  last  image  arrives  before  r0,  the 
diffraction  from  the  wedge  apex.  We  show  an  image  con¬ 
struction  on  Fig.  2.  The  images  are  reflections  of  the  source 
in  wedge  planes.  Possible  image  locations  are  on  the  circle 
containing  the  source  and  at  the  angles  given  by  2kSw±S0 
and  0=0.  The  choice  of  a  receiver  angle  0  requires  that  (8) 
be  tested  to  determine  the  active  images. 

The  rigid,  free,  and  free-rigid  boundaries  can  be  written 
as  follows 


/?(!)= -4tt[/,0  P0] 


Sit- [/?*(- -)/c]} 

+  2j  1,2  '™\<0  _ ) 


where  l  2  and  JBi  0  are  the  reflection  coefficients  at  the 
interfaces.  The  direct  arrival  is  given  by  k=0  and  the 
Rk(  -+)  term.  When  the  source  and  receiver  are  near  to  each 
other  or  at  the  same  place,  one  can  drop  this  term.  As  Chu 
has  shown,17,26  and  JBl  0  are  not  restricted  to  values  of 
±1,  but  can  have  any  values  between  —1  and  1  for  the 
isovelocity-density  contrast  wedge.  The  image  arrivals  are 
shown  in  Fig.  3. 

C.  Region  3:  Diffraction  and  t>r0 

The  final  closed  form  expressions  for  the  diffraction  ar¬ 
rivals  are  the  result  of  much  manipulation.  The  expressions 
use  the  following  parameters: 


77=  arc  cosh  U,  rf=ln[U+(U2-  1  )1/z], 

U=\c2t2  —  (r2  +  rl  +  y2)\/2r0r  ,  p=n/0w. 

Rigid  wedge  boundaries  (Biot-Tolstoy  and  Tolstoy)12,16,24 


P(0  =  — 4ir[P o  fl0] 


[Po  *o ]c 


[D(V,  +  -)+D(V,-  +  ) 


where 


20^0  1  '  7 
+  D(j7,  +  +)+D(  7, - )], 


sin  fr(7T+0-@o) 


D(v,+  )  sinh  j^cosh[ v( 77+ 5)]-cos  r'(ir+0-0o)} 


and  S  is  a  small  damping  number,  typically  10~5  or  a  small 
fraction  of  a  time  step. 

Free  wedge  boundaries  (Kinney  et  al.)25 


/?(/)= -4tt[P0  fl„] 


—  \.P o  R«]c 


[0(77,  +  -)+/>(  y,-  +  ) 


20wrro  L  v  '  v  n 
-0(77, +  +  )-/)( 77,--)]. 

Free-rigid  wedge  boundaries  (Chu)26 

d<f>  3 

p(t)=-4ir[P0  R0]-^- 
[P o  «o]c 

=  rroe°  \-Drlf(  77,  +  -  )  +  Or//(  77,  -  +  ) 
+  Of//(77,+  +)+Dr/f(Tj, - )]. 


Diffraction 

arrival 


1  2  3 


Time  (s) 

FIG.  3.  Signal  in  an  ideal  shallow  water  wedge  waveguide.  Signal  pressure 
p(t)  relative  to  the  amplitude  of  the  direct  arrival.  Source:  angle  8  deg  from 
surface  and  3000  m  range,  y  =0.  Receiver:  4  deg  from  surface  and  1000  m 
range,  y  =0.  The  input  signal  is  a  2-ms  spike. 
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Drlf(  V>+  ~  ) 

sin(  v!2){  ir+  0  -  @0)cosh[(  v/2)  y+  8/2] 
sinh  r^coshfH  *?+£)]“  cos  ^(7r  +  0-@o)}  ' 

A  sign  and  small  algebraic  error  in  Chu  Eq.  (14)  are 
corrected.26  For  special  values  of  0,, ,  the  diffraction  compo¬ 
nent  is  zero,  i.e., 

it  dd>\ 

if  ®iv=—  then  — —  =0,  (16) 

m  at 

where  m  is  an  integer.  The  images  are  the  complete  solution. 


D.  Numerical  evaluations  for  very  small  y 

These  expressions  (12)— (16)  are  easy  to  evaluate  for 
rf> 0  or  Tj  greater  than  the  first  time  step  beyond  r0  in  a 
numerical  evaluation.  As  t  goes  to  r0,  ry— *0,  and  dfajdt  be¬ 
comes  infinite.  One  can  average  D(A,+-)  over  the  small 
time  Tb  or  the  first  time  step  for  robust  evaluations. 

1  /*  7* 

D(A,+  -)|,  *D(A,  +  -)<fA,  A^r-r0, 

0  I  bJO 


where  A<Sr0.  Numerical  methods  are  in  references. 10,1 7,26,28 
Since  those  papers  were  written  we  have  become  ruthless  in 
dropping  small  terms,  because  only  the  first  time  step  is 
needed.  Simplified  results  follow. 

-  ^  2  sin  i/(ir+0-@o) 

,"‘T°  b^T^p^-ras^(rr+@^-0^Vr^/^  ’ 

r/^byjA,  b=c\l2r0/rr0.  (18) 

Similar  operations  on  the  free-rigid  wedge  (15)  give 

Dr/f(V,+  ~)\t~T0 

^  2  sin|>(?r+9-®oV2] 

~bjTh[l  —  cos  v(rr+0  — 0O)+ £2/2] 


An  example  of  numerical  evaluation  for  the  free-rigid  wedge 
is  shown  in  Fig.  3.  The  diffraction  arrival  is  the  last  event  in 
Fig.  3.  Comparisons  of  image  amplitudes  and  the  diffraction 
amplitude  shows  that  the  diffraction  carries  less  than  0.1%  of 
the  energy. 

Taking  a  practical  view  and  assuming  that  the  actual 
wedge  apex  is  a  lossy  beach,  the  wedge  apex  is  expected  to 
be  “leaky.”  Accordingly,  the  diffracted  pressure  component 
is  expected  to  be  negligible  compared  to  the  image 
contributions.28  We  ignore  the  diffraction  arrival  and  keep 
the  images. 


II.  IMAGING  OF  ACOUSTIC  FIELDS 

In  our  imaging  discussion  we  identify  the  transmission 
in  the  real  (actual)  wedge  as  having  the  subscript  A  and  the 
model  wedge  as  having  the  subscript  M.  Presumably  the 
model  wedge  is  an  accurate  model  of  the  actual  wedge.  We 
use  our  theoretical  solutions  to  compute  the  transmission 
from  a  source  any  place  in  the  model  wedge  to  a  probe 
receiver  at  position  Q. 


Briefly,  the  source  at  an  unknown  location  transmits  an 
impulsive  signal  and  the  transmission  is  received  by  a  hydro¬ 
phone  and  stored.  The  signal  is  time  reversed.  The  time  re¬ 
versed  signal  is  transmitted  into  the  model  wedge  by  a  trans¬ 
ducer  at  the  receiver  location.  Model  transmission  functions 
are  used  to  compute  the  pressures  at  a  set  of  locations.  A  map 
of  the  peak  pressures  gives  an  image  of  the  source  location. 

Recalling  (9)  and  (11),  the  model  transmission  functions 
can  be  defined  as  follows: 


fu(‘,  H,Q )- 


1*1  *(*-[**( -  +  )/c]} 
1,0  **(-+) 


|*-i| 

M  X  2  jB 


1*1  8{t-[Rk(--)/c]}\ 
1,0  **(--)  ) 


(20) 


and 


fM0,  H,Q )■  [D(V,  +  -)+D(V,-  +  ) 

+D(V,+  +  )  +  D(y,--)l  (21) 

where  the  subscript  H  refers  to  the  hydrophone  at  the  re¬ 
ceiver  location  (rH,  y //,©//)  and  Q  the  probe  location 
(rQ  >  '©q)  an(l  also  S  is  the  actual  source  location 

(rs,  ys,®s)-  Note  that  fM(t,  H,Q)  is  a  model  transmission 
function  and  does  not  include  the  source  function. 

In  the  actual  wedge,  the  pressure  signal  at  the  receiver 
for  a  transmission  of  s(t)  is  the  convolution 

ps(t)=s(t)*fA(t,S,H),  (22) 

where  *  indicates  the  convolution  operation  and  we  let  s(t ) 
be  the  impulsive  function  [P0  R0]<5(/).  The  time  reverse  is 

S,  H) 

=[P0  R0]S(-t)*fA{-t,S,H).  (23) 

Using  8(t)=  8(-t),  the  transmission  of  the  time  reversed 
signal  into  the  model  wedge  gives  the  pressure  at  the  probe 
position  Q, 

PmU,Q)  =  [P0  R0]S(t)*fA(-t,  S,  H)*fM(t,  H,Q). 

(24) 

The  image  is  formed  by  mapping  the  peak  of  | Pm((,Q)\  <n 
the  wedge  as  a  function  of  (rQ ,  yo,0e).  At  the  matched 
position 

cA(t,  S)  =  [P0  R0m*fA(-t,  S,  H)*fA(t,  H,  S). 

(25) 

The  operation  fA(~t,  S,  H)*fA(t,  H,  S)  gives  the  covari¬ 
ance  of  the  transmission  function.  We  call  this  the  correlation 
maximum. 

Figure  4(a)  shows  an  example  of  the  pressure  pM(t,  S ) 
at  the  matched  location  and  Fig.  4(b)  an  example  of  the 
pressure  Pu(t,Q)  at  an  unmatched  location.  In  these  figures, 
and  in  the  rest  of  the  paper,  the  source  and  receiver  angles 
and  'irH  etc.  are  measured  clockwise  from  the  top  face  of 
the  wedge  so  that 
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*S  =  ®tV-©5. 

'I'W=0H,-0H, 

and  (26) 

^0=©vv-e0. 

The  lateral  separation  y  is  measured  between  the  r-0  plane 
that  contains  the  source  and  the  r-0  plane  that  contains  the 
receiver.  An  experimental  verification  of  the  time  reverse 
transmission  technique  has  been  reported  in  the  well  known 
paper  by  Parvulescu  and  Clay/9 


III.  ENVIRONMENTAL  MODEL  AND  TRANSMISSIONS 

The  environmental  model  used  for  this  study  is  the 
three-dimensional  wedge  shown  schematically  in  Fig.  1.  We 
have  a  shallow  isospeed  water  layer  (c  =  1500  m/s)  retained 
within  the  wedge,  which  has  a  free  (pressure  release)  bound¬ 
ary  at  the  water  surface  and  a  rigid  bottom.  The  wedge  angle 
is  11  deg.  The  positions  of  the  source  and  receivers  within 
the  waveguide  are  shown  schematically  in  the  figure.  A  cy¬ 
lindrical  coordinate  system  is  used  to  specify  their  locations. 
The  variables  are  the  range  r  from  the  wedge  apex;  the  hori¬ 
zontal  distance  y  along  the  wedge  axis;  and  the  angle  of 
declination  'P  from  the  water  surface  boundary.  The  source  is 
placed  so  that  r5=3000  m,  ys=0  m,  and  lies  between  0 
and  11  deg.  The  receivers  are  placed  such  that  typically 
r/,^1000  m,  yH  takes  values  between  ±250  m,  and  lies 
between  0  and  11  deg.  One  of  the  features  of  this  study  is  the 
use  of  several  receivers  (displaced  from  each  other  by  vary¬ 
ing  rH,yH,  or  'Vw)  to  study  the  effect  on  localization  per¬ 
formance  of  multiple  receiver  channels.  If  only  one  receiver 
is  used  to  filter  the  signal  transmitted  by  the  source,  the  hy¬ 
drophone  is  positioned  at  rw=1000  m,  yH= 0  m,  and  'Pw=4 
deg. 

No  acoustic  attenuation  is  incorporated  into  the  propa¬ 
gation  paths  (except  1/r  spreading)  or  into  the  reflections  at 
the  wedge  boundaries.  Since  the  upper  wedge  boundary  is  a 
pressure  release  surface  (y8i  0=-l),  a  it  phase  change  is 
introduced  into  the  path  each  time  a  reflection  at  this  bound¬ 
ary  occurs.  This  causes  the  amplitude  of  the  signal  to  change 
sign.  In  Fig.  3  we  show  the  transmission  from  an  impulsive 
source  to  a  single  receiver.  The  input  signal  is  a  2-ms  boxcar 
approximation  to  a  delta  function  (a  spike).  The  source  is 
placed  at  rs=3000  m,  ys=0  m,  'Fs=8  deg.  The  first  two 
arrivals  are  seen  to  have  positive  amplitude.  These  are  due  to 
signals  transmitted  directly  from  the  source  and  from  the  first 
“  +  ”  image  of  the  source  reflected  in  the  bottom  surface  of 
the  wedge  where  «^li2=l  (see  Fig.  2).  The  next  four  arrivals 
have  negative  amplitude.  These  come  from  the  first  four 
“  -  ”  images  shown  in  Fig.  2.  All  of  the  ray  paths  here  in¬ 
volve  one  reflection  in  the  upper  surface  of  the  wedge  and 
therefore  undergo  a  v  phase  change.  There  are  then  another 
four  arrivals  with  positive  amplitude  (from  the  next  four 
“  +  ”  images).  These  ray  paths  undergo  two  reflections  in  the 
upper  surface  and  so  have  two  it  phase  changes  (i.e.,  a  total 
of  27t).  The  sequence  of  arrivals  continues  until  all  of  the 
images  shown  in  Fig.  2  are  accounted  for.  Inspection  of  Fig. 
3  shows  that  the  received  signal  has  about  1.3  s  duration.  The 
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FIG.  4.  Time  reversed  transmissions  to  matched  and  unmatched  locations. 
In  these  examples  the  receiver  location  is  4  deg  from  surface  and  1000  m 
range,  y~Q.  liie  matched  source  location  is:  angle  8  deg  from  surface  and 
3000  m  range,  y= 0.  The  unmatched  location  is  angle  8  deg  from  surface 
and  2750  m  range,  y=0.  The  input  signal  was  a  2-ms  spike,  (a)  Matched 
location,  (b)  Unmatched  location. 


backpropagation  is  calculated  by  time  reversing  the  received 
signal  and  retransmitting  it  from  the  receiver  (now  treated  as 
a  source),  and  from  all  the  images  of  the  receiver,  using  the 
same  formalism  as  for  the  forward  case.  An  example  of  the 
time-reversed  transmissions  to  a  matched  position  is  shown 
in  Fig.  4(a).  Spatial  resolution  in  source  location  depends  on 
the  reduction  of  the  correlation  maximum  as  the  receiver  is 
moved  away  from  the  matched  position.  Fig.  4(b)  shows  the 
effect  of  a  250-m  range  mismatch  for  the  transmission  shown 
in  Fig.  4(a). 

IV.  MEASURES  OF  LOCALIZATION  PERFORMANCE  ^ 

Localization  performance  depends  upon  the  number  and 
positioning  of  receivers  and  the  pulse  duration  (inversely 
proportional  to  the  signal  bandwidth).  Performance  is  repre¬ 
sented  by  the  relative  magnitudes  of  sideiobes  as  compared 
to  the  correlation  maximum  at  the  matched  position  (the 
main  peak)  and  the  spatial  resolution  of  the  peak. 

A.  Sideiobes 

A  maximum  locator,  defined  as  follows,  gives  a  simple 
measure  of  relative  amplitudes  of  the  sideiobes 
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FIG.  5.  Locator  performance  as  a  function  of  range  and  angle  of  declination.  In  this  example  the  locator  function  LP  defined  in  Eq.  (27)  is  plotted  for  a  single 
receiver  as  a  function  of  range  from  1500  to  4500  m,  and  declination  angle  from  I  deg  to  10  deg,  fory=0  m.  Source:  angle  8  deg  from  surface  and  3000  m 
range,  y  =0.  Receiver:  4  deg  from  surface  and  1000  m  range,  y  =0.  A  2-ms  input  spike  was  used. 


,  _max\pM(t,,Q)\ 

Lp~rmx\cA(tk,S)\' 

where  /,  and  tk  are  the  times  of  the  respective  maxima  of  the 
probe  signal  |pw(//,(?)|  and  the  auto  covariance  \cA(tk ,  S)|. 

A  plot  of  the  relative  maxima  of  Lp  as  the  probe  receiver 
is  moved  from  the  matched  position  is  a  simple  locator.  An 
image  of  Lp  for  a  single  receiver  is  shown  in  Fig.  5.  A  2-ms 
pulse  was  used  for  this  example.  The  source  was  again 
placed  at  rs=3000  m,  y  s=0  m,  = 8  deg.  The  value  of  Lp 
is  plotted  as  a  function  of  range  and  declination  angle.  The 
plot  shown  in  this  figure  clearly  indicates  a  main  peak  at  the 
correct  range  and  angle.  Spatial  resolution  is  very  good,  even 
though  only  a  single  receiver  has  been  used.  Surrounding  the 
main  peak  are  several  sidelobes.  These  are  more  pronounced 
on  the  side  closer  to  the  receiver.  Toward  the  apex  of  the 
wedge  they  become  quite  large  and  the  performance  of  the 
locator  deteriorates.  If  P  is  the  peak  level  in  Fig.  5,  while  p. 
is  the  mean  background  level  and  cr  the  standard  deviation  of 
the  background,  a  convenient  way  of  quantifying  the  peak- 
to-background  ratio  is  given  by 

PBR=10  1og10[(P-M)/ff]-  (28) 

This  can  be  viewed  as  a  measure  of  detectability  based  upon 
the  difference  between  the  peak  value  and  the  mean  back¬ 


ground  value,  scaled  in  units  of  the  standard  deviation  and 
expressed  in  dB. 

B.  Spatial  resolution 

How  may  the  spatial  resolution  of  the  source  in  Fig.  5  be 
quantified?  The  radial  width  of  the  main  signal  peak  is  a 
localization  parameter.  Figure  6(a)  shows  a  section  through 
the  peak,  taken  along  a  radial  from  the  wedge  apex,  for  the 
case  represented  in  Fig.  5.  We  see  that  the  maximum  of  Lp 
along  the  radial  has  a  form  very  similar  to  a  Gaussian,  with 
ripple  superimposed  on  it.  In  this  study  we  have  used  a  con¬ 
jugate  gradient  technique  to  fit  a  Gaussian  curve  to  this  data, 
so  that  the  signal  peak  radial  width  may  be  quantified  by  the 
Gaussian  half-width  p.  By  this  means  we  can  investigate 
how  the  addition  and  position  of  receivers,  and  changes  of 
pulse  duration,  broadens  or  sharpens  the  peak.  Another  im¬ 
portant  parameter  is  the  angular  width  of  the  main  signal 
peak.  Figure  6(b)  shows  a  section  through  the  signal  peak, 
taken  along  an  arc  centered  on  the  wedge  apex,  for  the  case 
represented  in  Fig.  5.  We  see  that  the  source  is  very  well 
resolved  in  angle,  even  with  only  one  receiver.  However,  a 
number  of  sidelobes  appear  at  discrete  angles  within  the 
wedge.  One  of  our  aims  here  to  demonstrate  how  these 
“angle”  sidelobes  may  be  eliminated  by  adding  more,  angu¬ 
larly  displaced,  receivers. 
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FIG.  6.  Locator  performance  along  radial  and  angular  sections  through  the 
peak.  Source:  angle  8  deg  from  surface  and  3000  m  range,  y=0.  Receiver: 
4  deg  from  surface  and  1000  m  range,  y  =0.  A  2-ms  input  spike  was  used, 
(a)  Lf  is  plotted  as  a  function  of  range  from  2750  to  3250  m,  while  the 
declination  angle  is  8  deg  and  y  =0.  (b)  Lf  is  plotted  as  a  function  of 
declination  angle  from  1  to  10  deg,  while  the  range  is  3000  m  and  y  =0. 


V.  SINGLE  RECEIVERS  AND  PULSE  DURATION 

The  main  peak  width  and  peak-to-background  ratios 
were  obtained  and  investigated  as  functions  of  (a)  single  re¬ 
ceiver  placement,  (b)  pulse  duration,  and  (c)  multiple  re¬ 
ceiver  placement.  For  reference  and  comparison  of  perfor¬ 
mance  of  single  and  multiple  receivers,  in  this  section  we 
describe  the  single  receiver  behavior. 

In  Fig.  7  we  plot  the  variation  of  PBR  and  p  for  a  2-ms 
pulse  case  in  which  a  source  is  placed  at  rs=3000  m,  y5=0 
m,  NPS=8  deg;  and  a  single  receiver  is  moved  such  that  rH 
varies  from  500-1500  m,  while  yH  and  are  maintained  at 
constant  values  of  0  m  and  4  deg,  respectively.  The  figure 
shows  that  both  quantities  improve  significantly  as  the  re¬ 
ceiver  is  moved  away  from  the  wedge  apex.  The  value  of 
PBR  increases  steadily  from  about  9  dB  when  rw=500  m  to 
about  11.6  dB  when  rH= 1250  m,  it  then  drops  slightly  to 
11.5  dB  at  rH= 1500  m.  The  Gaussian  half-width  p  falls  from 
240  m  to  about  45  m  over  the  same  interval.  The  reason  for 
this  improvement  in  performance  is  geometrical  and  easily 
understood  if  we  look  at  Fig.  2,  where  the  circle  of  source 
images  (these  will  be  images  of  the  receiver  in  the  reverse 
transmission  procedure)  is  constructed.  The  radius  of  this 
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FIG.  7.  Variation  of  peak-to-background  ratio  and  half-width  as  a  function 
of  single  receiver  placement.  PBR  and  p  are  plotted  for  a  source  at  angle  8 
deg  from  surface  and  3000  m  range,  y  =0.  A  2-ms  input  spike  was  used.  A 
single  receiver  was  placed  at  angle  4  deg  from  surface  and  y  =0.  The  re¬ 
ceiver  range  varied  from  500  to  1500  m. 


circle  increases  as  the  source  is  moved  further  from  the  apex. 
This  means  that  the  maximum  angle  subtended  by  tangents 
from  the  circle  to  any  point  within  the  wedge  also  increases, 
so  that  the  location  of  the  point  is  more  sharply  and  accu¬ 
rately  fixed  by  triangulation.  This  decreases  p  directly,  and 
also  reduces  the  background  variance  and  increases  PBR  by 
decreasing  the  opportunity  for  constructive  addition  of  sig¬ 
nals  on  different  ray  paths  at  locations  other  than  the 
matched  position. 

The  angular  placement  of  the  receiver  can  also  have  a 
very  noticeable  effect  on  the  performance  of  the  locator.  In 
Fig.  8(a)  we  show  the  variation  of  Lp  along  an  arc  at 
r (2 =1000  m  for  a  2-ms  pulse  case  in  which  the  source  is 
placed  at  rs=3000  m,  ys=0  m,  'P5=5.5  deg;  and  the  re¬ 
ceiver  is  placed  at  rH=  1000  m,  yH= 0  m,  ,RW=4  deg.  The 
main  signal  peak  is  seen  at  the  correct  source  angle.  How¬ 
ever,  there  are  also  two  prominent  sidelobes  positioned  sym¬ 
metrically  on  either  side  of  the  main  peak.  In  Fig.  8(b)  we 
show  the  equivalent  curve  when  the  receiver  angle  is 
tyH=8.5  deg,  while  all  the  other  parameters  are  unchanged. 
Again  the  main  peak  is  clearly  distinguished,  but  the  side- 
lobes  have  disappeared.  The  sidelobe  angles  in  Fig.  8(a) 
seem  to  be  related  to  the  relative  angles  of  the  source  and 
receiver.  If  we  write  ^s=5.5  deg  and  'Pw=4  deg,  then  it 
appears  =  'PS±2{|,P5-'P//|}  (giving  VL=2.5  and  8.5 
deg  in  this  case)  for  the  sidelobes.  When  ^w=8.5  deg  the 
same  formula  gives  'Pi  =  -0.5  and  11.5  deg,  placing  the 
sidelobes  outside  of  the  wedge  so  that  they  do  not  appear  in 
Fig.  8(b). 

The  effect  of  bandwidth  on  the  performance  of  the  loca¬ 
tor  may  be  investigated  indirectly  by  varying  the  time  dura¬ 
tion  of  the  pulse  signal.  In  the  cases  we  have  looked  at  so  far 
we  have  used  a  2-ms  boxcar  spike.  This  corresponds  to  a 
bandwidth  si  kHz.  Now  we  will  introduce  a  1-ms  spike 
(s2-kHz  bandwidth)  and  a  4-ms  spike  (=500- Hz  band- 
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FIG.  8.  Locator  performance  along  angular  sections  for  different  receiver 
angles.  The  source  is  placed  at  angle  5.5  deg  from  surface  and  3000  m 
range,  y  =0.  A  2-ms  input  spike  was  used.  LF  is  plotted  as  a  function  of 
declination  angle  from  1  to  10  deg,  while  the  range  is  3000  m  and  y  =0.  (a) 
The  receiver  is  placed  4  deg  from  surface  and  1000  m  range,  y=0.  (b)  The 
receiver  is  placed  8.5  deg  from  surface  and  1000  m  range,  y  =0. 

width)  also.  In  Fig.  9(a)  we  show  the  variation  of  p  for  the 
three  pulse  lengths,  for  a  case  in  which  a  source  is  placed  at 
rs= 3000  m,  ys=0  m,  ^=8  deg;  and  a  single  receiver  is 
moved  such  that  rH  varies  from  500-1500  m,  while  yH  and 

are  maintained  at  constant  values  of  0  m  and  4  deg, 
respectively  (this  is  the  same  deployment  used  to  produce 
Fig.  7).  The  figure  clearly  shows  that  the  spatial  resolution  of 
the  main  peak,  as  measured  by  p,  improves  in  direct  propor¬ 
tion  to  the  bandwidth  and  inversely  to  the  duration  of  the 
pulse.  The  values  of  p  shown  in  the  figure  do  not  indicate  a 
strict  linear  relationship  with  the  pulse  duration.  Figure  9(a) 
shows  the  narrowing  behavior  of  the  main  peak  for  all  three 
pulses  as  the  receiver  is  moved  further  from  the  apex.  In  Fig. 
9(b)  we  plot  the  PBR  values  for  the  three  pulses  as  a  function 
of  receiver  location.  This  quantity  again  improves  in  propor¬ 
tion  to  the  bandwidth  and  clearly  increases  as  the  receiver  is 
moved  further  from  the  apex,  although  some  small  fluctua¬ 
tions  may  be  observed  for  the  1-  and  2-ms  pulses.  Doubling 
the  bandwidth  seems  to  give  between  0.25-1.5  dB  of  en¬ 
hancement  in  PBR  at  each  receiver  location. 

VI.  MULTIPLE  RECEIVER  ARRAY 

The  rationale  for  putting  multiple  receiver  arrays  into 
the  waveguide  is  that  additional  information  is  made  avail- 
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FIG.  9.  Variation  of  peak-to-background  ratio  and  half-width  as  a  function 
of  pulse  duration.  The  source  is  placed  at  angle  8  deg  from  surface  and  3000 
m  range,  y  =0.  A  single  receiver  was  placed  at  angle  4  deg  from  surface  and 
y=0.  The  receiver  range  varied  from  500  to  1500  m.  Three  pulses  were 
used:  1,  2,  4  ms.  (a)  Peak  half-width  p.  (b)  PBR. 

able  about  the  source  by  virtue  of  the  multiplicity  of  ray 
paths  to  the  many  receiver/image  locations.  Arrays  are  ex¬ 
pected  to  improve  spatial  resolution,  reduce  the  background 
variance,  and  remove  ambiguities  by  enhancing  the  effec¬ 
tiveness  of  source  localization  by  triangulation,  as  we  ob¬ 
served  earlier  while  studying  the  placement  of  a  single  re¬ 
ceiver.  When  adding  further  receivers,  we  can  vary  their 
placement  position  by  changing  rH,  yH,  or  'PH .  These  are 
orthogonal  parameters  within  the  wedge,  and  we  will  inves¬ 
tigate  their  effect  on  the  localization  process  independently 
of  each  other.  Just  as  we  did  earlier  in  the  study  of  single 
receiver  position  we  will  quantify  this  process  by  measuring 
the  values  of  PBR  and  p,  this  time  as  functions  of  the  array 
length. 

In  Fig.  10(a)  we  plot  the  variation  of  PBR  and  p  for  a 
2-ms  pulse  where  the  source  is  placed  at  rs=3000  m,  _y5=0 
m,  ¥$=8  deg.  Initially  a  single  receiver  is  introduced  at 
rw=1000  m,  yH— 0  m,  ^=4  deg.  Additional  receivers  are 
positioned  so  that  they  lie  along  a  radial  passing  through  the 
first  receiver.  The  second  is  placed  at  rH= 950  m,  yH= 0  m, 
,'PW=4  deg,  forming  a  two  element  array  of  length  50  m.  The 
third  is  placed  at  rw=1050  m,  yH= 0  m,  ^=4  deg,  making 
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FIG.  10.  Variation  of  peak-to-background  ratio  and  half- width  as  a  function 
of  array  length.  The  source  is  placed  at  angle  8  deg  from  surface  and  3000 
m  range,  y  =0,  A  2-ms  input  spike  was  used.  The  first  receiver  is  placed  at 
angle  4  deg  from  surface,  range  1000  m  and  y  =0.  (a)  PBR  and  p  are  plotted 
as  additional  receivers  are  introduced  at  SO  m  range  intervals,  with  the  same 
angle  and  y  as  the  first  receiver,  (b)  PBR  and  p  are  plotted  as  additional 
receivers  arc  introduced  at  1-deg  angular  intervals,  with  the  same  range  and 
y  as  the  first  receiver. 


a  100-m  array.  A  fourth  and  fifth  receiver  are  introduced  at 
rw=900  m  and  rH=  1100  m,  increasing  the  array  length  to 
150  and  200  m,  respectively.  We  see  that  with  only  one  re¬ 
ceiver  the  Gaussian  half-width  p  is  about  67.5  m.  When  the 
receiver  at  950  m  is  added,  p  increases  to  slightly  less  than 
75  m.  The  reason  for  this  increase  is  clear  from  Fig.  7.  Re¬ 
ceivers  closer  to  the  wedge  apex  resolve  the  source  less  well, 
so  that  backpropagation  from  the  second  receiver  produces  a 
broader  signal  peak  than  from  the  first.  Locating  on  the  ag¬ 
gregate  signal  from  the  two  receivers  (i.e.  processing  them  as 
a  two  element  array)  gives  a  poorer  result  than  is  achieved 
with  the  first  receiver  alone.  The  opposite  effect  is  observed 
when  a  third  receiver  is  added  at  1050  m,  which  decreases  p 
to  68  m.  The  pattern  is  repeated  when  the  fourth  and  the  fifth 
receivers  are  introduced:  p  increases  to  71  m  and  then  de¬ 
creases  again  to  66  m.  It  appears  that,  contrary  to  expecta¬ 
tions,  the  aggregation  of  signals  from  multiple  receivers  does 
not  lead  to  a  monotonic  improvement  in  spatial  resolution  of 
the  source,  unless  the  secondary  receivers  are  further  from 


FIG.  11.  Suppression  of  "angle”  sidelobes  by  addition  of  receivers.  The 
source  is  placed  at  angle  5.5  deg  from  surface  and  3000  m  range,  y  =0.  A 
2-ms  input  spike  was  used.  All  receivers  are  placed  at  range  1000  m  and 
y  =0.  LP  is  plotted  as  a  function  of  declination  angle  from  1  to  10  deg,  while 
the  range  is  3000  m  and  y  =0.  (a)  Two  receivers:  3  and  4  deg.  (b)  Three 
receivers:  3,  4,  and  5  deg.  (c)  Four  receivers:  2,  3,  4,  and  5  deg. 

the  wedge  apex  than  the  primary  receiver.  This  is  a  surpris¬ 
ing  result,  but  more  significant  is  the  fact  that  the  fluctuations 
in  p  aiv  so  small.  The  spread  in  values  is  only  9  m;  and  the 
200-m  long  array  has  the  total  effect  of  reducing  p  by  only 
1.5  m  from  the  value  achieved  using  a  single  receiver.  This  is 
in  complete  contrast  to  the  improvements  that  can  be  ob¬ 
tained  by  judicious  positioning  of  a  single  receiver,  as  we 
observed  earlier.  The  variations  in  PBR  shown  in  the  same 
figure  mirror  and  complement  the  p  results.  This  falls  from 
an  initial  value  of  10.9  dB  for  the  single  (1000  m)  receiver  to 
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FIG.  12.  Locator  performance  as  a  function  of  range  and  y  for  multiple  lateral  receivers.  The  locator  function  L P  is  plotted  as  a  function  of  range  from  2000 
to  4000  m,  and  y  from  - 1000  to  1000  m,  with  a  declination  angle  of  4  deg.  Source:  angle  8  deg  from  surface  and  3000  m  range,  y=350  m.  A  2-ms  input 
spike  was  used.  All  receivers  are  placed  4  deg  from  the  surface  and  at  range  1000  m.  (a)  One  receiver:  >’  =  150  tn.  (b)  Two  receivers:  >>  =  150  and  200  m. 
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10.8  dB  when  the  second  (950  m)  receiver  is  added.  It  main¬ 
tains  this  value  with  the  third  (1050  m)  receiver  and  then 
falls  again  to  10.6  dB  with  the  fourth  (900  m).  Finally  it  rises 
slightly  to  10.7  dB  with  the  introduction  of  the  fifth  (1100  m) 
receiver.  All  of  the  variations  are  very  small,  with  a  spread  of 
only  0.3  dB.  It  seems  that  the  addition  of  more  receivers,  and 
aggregating  signals  over  a  wide  range  aperture,  does  not 
smooth-out  sidelobes  in  the  waveguide  or  decrease  the  back¬ 
ground  variance.  Again,  it  would  be  better  to  use  a  single 
receiver  and  move  it  further  from  the  wedge  apex.  We  also 
looked  at  the  change  in  angular  resolution  of  the  source  with 
the  use  of  multiple,  radially  displaced,  receivers.  No  signifi¬ 
cant  effect  was  observed. 

In  Fig.  10(b)  we  again  plot  the  variation  of  PBR  and  p 
for  a  2-ms  pulse  with  the  source  placed  at  r5=3000  m,  ys=0 
m,  'Pv = 8  deg.  As  before,  the  first  receiver  is  placed  at 
rw  =  1000  m,  yH  =  0  m,  'irH=4  deg.  This  time,  however,  ad¬ 
ditional  receivers  are  positioned  so  that  they  lie  along  an  arc 
passing  through  the  first  receiver.  The  second  is  placed  at 
rw=1000  m,  yH= 0  m,  ^=3  deg,  forming  a  two  element 
array  of  angular  aperture  1  deg.  The  third  is  placed  at 
rH  =  1000  m,  yw=0  m,  'Fw=5  deg,  making  a  2-deg  array.  A 
fourth  and  fifth  receiver  are  introduced  at  'Pw=2  deg  and 
'PH= 6  deg,  increasing  the  array  aperture  to  3  and  4  deg, 
respectively.  With  only  one  receiver  the  Gaussian  half-width 
p  is  about  67.5  m.  When  the  receiver  at  ^H=3  deg  is  added, 
p  increases  to  slightly  less  than  74  m.  It  then  decreases  again 
as  the  'Pw=5  deg  is  added  and  levels  off  at  around  70  m.  As 
in  Fig.  10(a),  the  PBR  curve  appears  to  mirror  the  p  curve.  It 
falls  from  an  initial  value  of  10.9  dB  to  about  10.3  dB  for 
two  receivers.  It  then  increases  again  to  a  peak  of  11.1  dB, 
and  falls  finally  to  10.9  dB  for  the  4-deg  aperture.  These 
changes  in  p  and  PBR  are  probably  due  to  the  introduction  of 
receivers  which,  depending  upon  the  contingencies  of  their 
own  location,  give  better  or  poorer  individual  performance. 
Again  they  are  quite  minor,  and  it  seems  that  the  addition  of 
more  receivers,  and  aggregating  signals  over  the  angular  ap¬ 
erture,  does  not  smooth-out  sidelobes  in  the  waveguide  or 
decrease  the  background  variance. 

Although  the  introduction  of  angularly  displaced  receiv¬ 
ers  does  not  change  the  radial  width  p  of  the  peak,  it  has  a 
notable  effect  on  the  angular  resolution  performance.  In  Fig. 
11(a)  we  show  the  variation  of  Lp  along  an  arc  at  rG  =  1000 
m  for  a  2-ms  pulse  with  the  source  placed  at  rs=3000  m, 
y5=0  m,  ^5=5.5  deg.  Two  receivers  are  placed:  one  at 
rw=1000  m,  yH= 0  m,  'Pw=4  deg;  and  a  second  with  ¥w= 3 
deg.  This  figure  should  be  compared  with  Fig.  8(a),  which 
was  produced  using  the  same  deployment  but  with  only  the 
4^=4  deg  receiver.  We  see  that  the  angular  resolution  of  the 
main  peak  is  about  the  same  in  the  two  cases,  but  the  ampli¬ 
tude  of  the  sidelobes  on  either  side  of  the  peak  is  much 
reduced  due  to  the  presence  of  the  second  receiver.  In  Figs. 
11(b)  and  11(c)  we  show  the  response  when  receivers  are 
added  at  ¥ H=5  deg  and  2  deg,  respectively.  In  Fig. 
11(b)  the  sidelobes  may  still  be  discerned,  but  in  Fig.  11(c) 
they  have  essentially  disappeared. 

Let  us  now  look  at  the  effect  of  positioning  multiple 
receivers  at  laterally  displaced  locations  (i.e.  with  variable 
yH  values).  In  Fig.  12(a)  we  plot  Lp  as  a  function  of  rQ  and 


FIG.  13.  Variation  of  peak-to-background  ratio  with  number  of  lateral  re¬ 
ceivers.  Source:  angle  8  deg  from  surface  and  3000  m  range,  y= 350  m.  A 
2-ms  input  spike  was  used.  All  receivers  are  placed  4  deg  from  the  surface 
and  at  range  1000  m.  PBR  is  plotted  as  receivers  are  introduced  at  y  =  150, 
200,  100,  250,  and  50  m. 

yQ  (with  Wq  held  constant  at  4  deg)  for  a  2-ms  pulse  case 
with  the  source  at  rs=3000  m,  _ys=350  m,  '4's=8  deg  and  a 
single  receiver  at  r//  =  1000  m,  yw=150  m,  4^=4  deg.  We 
see  that  the  one  receiver  is  able  to  locate  the  source  quite 
sharply,  but  that  it  gives  rise  to  mirror-image  ambiguities  on 
either  side  of  a  radial  passing  through  the  receiver  location  at 
yG  =  150  m,  so  that  a  second  signal  peak  appears  atyG  =  — 50 
m.  The  introduction  of  a  second  receiver  at  yw  =  200  m  im¬ 
mediately  removes  the  ambiguity  and  improves  the  spatial 
resolution,  as  we  see  in  Fig.  12(b).  The  addition  of  more 
receivers  also  reduces  the  background  variance.  Figure  13 
shows  that  PBR  increases  with  the  number  of  receivers,  in¬ 
dicating  enhanced  smoothing  of  sidelobes,  although  the  rate 
of  improvement  seems  to  slowly  dec.  ease.  We  also  studied 
the  effect  of  adding  laterally  displaced  receivers  on  the  radial 
width  (p)  and  the  angular  resolution  of  the  source  peak.  No 
significant  effect  was  observed. 

VII.  CONCLUSIONS 

Adding  more  receivers,  in  either  a  radial  or  angular  con¬ 
figuration,  has  an  insignificant  effect  in  reducing  the  peak 
width  or  improving  the  peak-to-background  ratio.  The  addi¬ 
tion  of  more  receivers  in  a  lateral  configuration  leads  to  sharp 
and  unambiguous  azimuthal  localization  of  the  source  and 
improves  the  peak-to-background  ratio.  The  range-dependent 
characteristics  of  a  wedge  waveguide  facilitate  source  local¬ 
ization.  “Range”  sidelobes,  which  are  a  typical  complication 
in  a  range-independent  environment,  are  eliminated  even 
when  only  a  single  hydrophone  is  used.  “Angle”  sidelobes 
may  be  effectively  eliminated  by  introducing  additional  re¬ 
ceivers  angularly  spaced  in  the  waveguide,  at  a  constant 
range  from  the  apex.  Both  the  peak  width  and  the  peak-to- 
background  ratio  improve  in  proportion  to  the  inverse  band¬ 
width  of  the  input  signal.  If  one  has  the  choice  to  make 
location  experiments,  wedges  are  much  better  than  nice  uni¬ 
form  waveguides.  It  may  be  speculation,  but  we  think  that 
source  localization  improves  with  the  complexity  of  the 
range  dependency  of  the  environment.  Of  course,  this  re¬ 
quires  solution  of  the  acoustic  transmission  problem  for  a 
large  frequency  bandwidth. 
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APPENDIX:  GEOMETRICAL  CONSTRUCTION  OF 
REFLECTED  RAY  PATHS 

The  total  pressure  field  at  the  receiver  location  due  to  the 
direct  arrival  from  the  source  and  the  reflected  arrivals  from 
the  images  is  given  by  Eq.  (9),  which  contains  two  summa¬ 
tions.  Examination  of  th<s  equation  shows  that  each  separate 
term  within  the  summations  uniquely  represents  a  single 
transmission  path  from  either  the  source  or  one  of  its  images 
to  the  receiver;  and  that  each  term  also  incorporates  (via  >0 
and  the  boundary  reflections  experienced  by  the  trans¬ 
mission  as  it  traverses  the  path  which  it  represents.  We  also 
note  that  the  source  and  images  are  represented  by  terms  in 
either  the  (-+)  summation  or  the  ( — )  summation,  indicat¬ 
ing  that  the  transmissions  fall  into  two  distinct  families. 

A  simple  and  beautiful  geometrical  construction  may  be 
performed  which  clearly  illustrates  the  application  of  Eq.  (9) 
and  facilitates  the  calculation  of  the  ray  paths,  especially  in 
cases  where  the  reflection  coefficients  J%10  and  .'Ml  2  are  de¬ 
pendent  upon  the  angle  of  incidence  at  the  boundaries.  It  also 
greatly  simplifies  transmission  calculations  in  a  finite  wedge, 
as  we  shall  also  show.  To  demonstrate  this  construction  we 
consider  the  case  of  a  35-deg  wedge,  with  a  source  and  re¬ 
ceiver  positioned  so  that  the  source  is  closer  to  the  apex. 
Figure  A1  shows  the  wedge  and  the  source  and  receiver  lo¬ 
cations,  and  uses  several  images  of  the  source  to  construct 
one  of  the  ten  or  so  available  ray  paths  from  the  source  to  the 
receiver.  We  see  in  the  figure  the  source  (labeled  A)  and  all 
the  images  of  the  source  symmetrically  positioned  around 
the  apex  as  determined  by  the  image  reconstruction  theory 
described  in  Sec.  I.  The  particular  ray  depicted  originates  at 
A  and  is  first  reflected  by  the  bottom  at  B,  making  an  angle 
aB  with  the  bottom  surface.  (For  the  purposes  of  this  discus¬ 
sion  it  is  convenient  to  use  the  complement  to  the  angle  of 
incidence  rather  than  the  angle  of  incidence  itself.)  When 
reflected  the  ray  is  modified  by  the  reflection  coefficient 
which,  if  it  is  angle  dependent,  may  be  written  in  this  case  as 
JB i_2  («b)-  The  ray  then  travels  to  the  water  surface  which  it 
meets  at  C.  It  makes  an  angle  Oq  with  this  boundary  and 
reflects  from  it,  modifying  the  ray  by  (oq).  The  ray  then 
reflects  a  second  time  from  the  bottom  surface  at  D:  the 
angle  is  aD  and  the  reflection  coefficient  is  J8lf2  (aD).  Finally 
the  ray  reaches  the  receiver  at  F.  The  partial  pressure  field  at 
F  due  to  this  single  ray  may  be  determined  by  isolating  the 
appropriate  term  from  the  summations  in  Eq.  (9).  In  this  case 
it  may  be  written 

PA—ri{)  =  [^*0  *o]^l ,2( ,o( ,2( ®d) 


^f-R_,(--)/c] 

*-i(— ) 


(Al) 


k  =  2 


FIG.  Al.  A  sample  ray  path  construction.  Several  of  the  source  images  are 
used  to  construct  one  of  the  possible  ray  paths  between  the  acoustic  source 
at  A  and  the  receiver  at  F.  A  similar  construction  may  be  performed  for  each 
of  the  images  shown.  The  images  of  the  surfaces  of  the  wedge  reflected  in 
each  other  are  also  shown.  The  “kaleidoscopic”  relationship  between  the 
source  and  surface  images  is  apparent.  Geometrical  reasoning  may  be  used 
to  demonstrate  the  equivalence  of  all  the  angles  labeled  aB ,  and  also  those 
labeled  ac  and  aD . 


We  see  that  k= — 1  for  this  transmission  and  that  it  lies  in  the 
( — )  family.  The  total  length  of  the  ray  path  is  /? _j( — ). 

In  addition  to  the  images  of  the  source  reflected  in  the 
surfaces  of  the  wedge.  Fig.  Al  also  shows  the  images  of  the 
wedge  surfaces  themselves  reflected  in  each  other.  (Careful 
examination  of  the  figure  shows  that  the  surfaces  are  re¬ 
flected  and  rereflected  in  each  other  to  form  a  kaleidoscope 
of  images  around  the  apex.)  We  may  describe  the  pattern  that 
results  from  this  as  follows:  Starting  at  the  water  surface  and 
proceeding  in  a  counterclockwise  direction  around  the  apex 
we  encounter,  after  turning  35  deg,  the  first  (k=l)  image  of 
the  bottom  surface  reflected  in  the  water  surface  (this  is 
drawn  as  a  thick  line  with  long  dashes).  Turning  another  35 
deg  we  find  the  first  image  of  the  water  surface  (drawn  as  a 
thinner  line  with  shorter  dashes).  Another  35  deg  brings  the 
second  {k=2)  image  of  the  bottom  surface,  and  a  further  35 
deg  leads  to  the  second  image  of  the  water  surface.  If  we 
then  go  over  to  the  original  bottom  surface,  and  start  turning 
successive  increments  of  35  deg  in  a  clockwise  direction,  we 
can  explain  the  remaining  surface  images  in  a  similar  man¬ 
ner.  We  should  note  here  the  angular  relationship  of  these 
wedge  surface  images  to  each  other  and  to  the  images  of  the 
source.  Bearing  in  mind  the  mathematical  constructions  of 
both  types  it  follows,  after  some  consideration  of  the  reflec¬ 
tion  symmetry  in  the  problem,  that  the  images  of  the  surfaces 
perfectly  bisect  the  angular  locations  of  the  pairs  of  source 
images  between  which  they  pass.  (It  again  helps  to  think  in 
terms  of  a  kaleidoscope  to  see  this.)  Using  basic  geometrical 
arguments  the  reader  may  easily  satisfy  himself  that  this  is 
true  and  also  prove  the  equivalence  of  the  sets  of  angles 
labeled  aB  and  oq  in  Fig.  Al. 

In  Fig.  A2  we  consider  the  same  ray  transmission  de¬ 
picted  in  Fig.  Al,  but  this  time  it  is  represented  by  a  straight 
line  drawn  directly  between  the  originative  source  image  (la¬ 
beled  A')  and  the  receiver  at  F.  This  is  an  equivalent  con¬ 
struction  of  the  ray  path.  We  can  show  this  and  also  demon- 
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FIG.  A2.  Alternative  ray  path  construction.  An  equivalent  description  of  the 
ray  path  in  Fig.  A1  is  shown.  Reflection  at  the  surfaces  of  the  wedge  is  now 
represented  by  transmission  through  the  images  of  those  surfaces  with  the 
same  angles  of  incidence  (compare  angles  trB ,  etc ,  and  aD  in  this  figure  and 
Fig.  Al).  The  total  path  length  and  the  lengths  of  the  various  sections  of  the 
paths  are  the  same  in  both  figures.  The  ray  path  expansion  [Eq.  (Al)]  con¬ 
tains  the  product  of  three  reflection  coefficients.  The  image  originating  the 
ray  is  in  the  ( — )  family. 


strate  that  all  of  the  events  occurring  within  the  transmission 
path  from  A  to  F  shown  in  Fig.  Al  may  be  correctly  (and 
more  easily)  described  by  considering  the  transmission  path 
from  A'  to  F  in  Fig.  A2.  First  of  all  we  will  consider  the 
overall  length  of  the  two  paths.  By  comparing  Figs.  Al  and 
A2,  and  using  symmetry  and  geometry,  we  can  see  that  the 
line  segment  AB  is  the  same  length  as  A'B'.  Similarly  BC  is 
the  same  length  as  B'C',  and  CD  is  the  same  length  as  C'D'. 
The  segments  DF  and  D'F  are  identical.  Adding  the  seg¬ 
ments  together  for  the  two  figures  shows  that  the  two  paths 
are  the  same  length  and  equal  to  f?_,( — ).  Now  let  us  fol¬ 
low  the  transmission  from  A'  to  F  in  Fig.  A2.  The  ray  starts 
at  A'  and  travels  a  distance  A'B'  (same  as  AB)  to  the  k = - 1 
image  of  the  bottom  surface,  which  it  meets  with  angle  aB . 
Since  this  surface  image  is  a  reflection  of  the  ^12  boundary, 
passage  of  the  ray  through  it  indicates  a  modification  of  the 
ray  by  the  reflection  coefficient  J8l  2(aB).  The  ray  then  trav¬ 
els  B'C'  (same  as  BC)  to  the  next  image  of  the  water  sur¬ 
face,  which  it  meets  with  angle  ate .  This  image  is  a  reflection 
of  the  0  boundary,  and  transmission  through  it  therefore 
indicates  a  modification  of  the  ray  by  j0(Qt)-  The  ray  then 
travels  C'D'  (same  as  CD)  to  the  bottom  surface,  which  it 
meets  with  angle  aD .  Passage  through  this  indicates  a  modi¬ 
fication  of  the  ray  by  J^12(aD).  The  ray  then  travels  D'F  to 
the  receiver.  If  we  make  a  piece  by  piece  comparison  of  the 
path  just  described  with  the  path  AF  in  Fig.  Al  we  see  that, 
at  every  stage,  they  are  identical  and  are  both  described  by 
Eq.  (Al).  The  partial  pressure  at  the  receiver  F  due  to  the  ray 
A'F  is  given  by  pA'->f(0=Fa-f(0-  The  maul  difference  in 
the  usage  of  Eq.  (Al),  when  working  in  the  image  space 
shown  in  Fig.  A2,  is  that  the  changes  to  the  ray  indicated  by 
the  reflection  coefficients  ,  J*?1>2  etc.  are  applied  as  the 
ray  passes  through  the  surfaces,  rather  than  being  reflected 
by  them. 

Using  Fig.  A 2  it  is  possible  to  provide  further  illumina¬ 


FIG.  A3.  A  higher-order  ray  path.  The  ray  shown  here  involves  fotr  reflec¬ 
tions  at  the  wedge  surfaces,  and  this  is  represented  by  transmission  through 
the  surface  images  as  indicated.  The  ray  path  expansion  [Eq.  (A2)]  contains 
the  product  of  four  reflection  coefficients.  The  image  originating  the  ray  is  in 
the  (-+)  family. 


tion  of  Eq.  (Al)  and  of  all  the  other  terms  in  Eq.  (10).  We 
recall  that  the  ray  path  AF  in  Fig.  Al  (and  hence  A'F  in  Fig. 
A2)  lies  in  the  ( — )  family  of  terms  with  k  =  —  1.  Looking  at 
Fig.  A2  we  note  that  the  source  image  A'  is  positioned  im¬ 
mediately  beside  the  it = - 1  reflection  of  the  bottom  surface. 
It  turns  out  that  the  ray  path  to  the  receiver  from  the  image 
A",  on  the  opposite  side  of  the  k=  -1  bottom  surface  reflec¬ 
tion  from  A',  is  also  given  by  a  k=  - 1  term  in  Eq.  (10),  but 
lies  in  the  (-+)  family.  This  is  true  for  all  the  pairs  of 
images  lying  on  either  side  of  the  various  reflections  of  the 
bottom  surface  in  the  image  space.  The  ray  paths  from  them 
to  the  receiver  will  be  labeled  by  the  k  value  of  the  bottom 
reflection  to  which  they  are  adjacent,  and  they  will  alter¬ 
nately  fall  into  the  ( — )  and  (- +)  families.  We  can  now  see 
that,  rather  than  thinking  in  terms  of  the  ray  paths,  we  can 
identify  the  images  themselves  according  to  their  k  value  and 
family.  In  Fig.  A3  we  show  a  ray  path  between  another  of  the 
source  images  and  the  receiver.  This  image  lies  adjacent  to 
the  k=2  reflection  of  the  bottom  surface  and  is  identified  as 
a  k—2  image  in  the  (-  +)  family.  A  straight  line  drawn  be¬ 
tween  it  and  the  receiver  indicates  that  the  ray  undergoes 
four  boundary  interactions  this  time  (note  that  in  this  figure 
we  have  used  a  fresh  lettering  scheme  for  the  source  and 
receiver  and  points  of  interaction  with  the  boundaries).  The 
partial  pressure  field  may  again  be  determined  by  isolating 
the  appropriate  term  from  the  summations  in  Eq.  (10).  In  this 
case  it  is 

P A— f( 0-[Po  ^o]^1  ,2( “b)*^! ,o( arc)-^l,2( aD)-^l ,o( <*e) 


S{t-R2(-  +)/c] 
*2(-+) 


(A2) 


From  the  foregoing  discussion,  we  can  now  see  how 
algorithms  may  be  developed  to  perform  ray  transmission 
computations  for  the  general  wedge  problem,  in  which  the 
wedge  angle  as  well  as  the  source  and  receiver  locations  and 
the  wedge  surface  reflectivities  are  variable  parameters.  The 
ray  paths  may  be  calculated  geometrically.  First,  the  wedge 
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FIG.  A4.  Allowed  and  disallowed  ray  paths  in  a  finite  wedge.  Images  of  the 
source  and  surfaces  are  shown  in  a  wedge  with  finite  faces  of  unequal 
length.  The  ray  path  AF  will  successfully  reflect  between  the  two  surfaces  of 
the  wedge  and  connect  the  source  and  receiver.  The  other  path  A'F  misses 
one  surface  at  D'  and  therefore  does  not  exist. 

would  be  drawn  and  all  the  reflections  of  the  bottom  and 
water  surfaces.  Second,  the  source  would  be  positioned  and 
all  of  its  images  together  with  the  receiver.  Third,  rays  would 
be  drawn  between  the  images  and  the  receiver,  and  the  path 
lengths  Rk( — )  and  Rk(—+)  calculated  for  all  k.  Fourth,  the 
angles  at  which  the  rays  cross  the  surface  reflections  ( aB , 
ac,...,  etc)  and  the  associated  reflectivities  [.^12(«b)’ 
0(ac),...  etc]  would  be  determined  (we  can  represent  the 
rays  and  surfaces  algebraically,  so  finding  their  points  and 
angles  of  intersection  is  straightforward).  Finally,  the  pres¬ 
sure  field  contributions  from  all  of  the  images  may  be  enu¬ 
merated  individually  [Eqs.  (Al)  and  (A2)  give  two  ex¬ 
amples]  and  then  added  together. 

A  particularly  interesting  consequence  of  this  procedure 
applies  to  the  case  of  finite  wedges,  i.e.,  when  the  faces  are 
of  finite  and  generally  unequal  length.  In  Fig.  A4  we  show  a 
35-deg  wedge  for  which  the  lower  surface  is  much  shorter 
than  the  upper  surface.  An  important  question  to  ask  here  is, 
for  an  arbitrary  placement  of  source  and  receiver,  which  ray 
paths  will  be  able  to  perform  multiple  reflections  from  both 
surfaces  of  the  wedge  and  successfully  connect  the  source 
and  receiver?  This  question  is  quite  simply  answered  if  we 
construct  the  images  of  the  source  around  the  apex,  as  be¬ 
fore,  and  also  the  images  of  the  surfaces.  We  note  that  the 
surface  images  are  now  of  different  lengths;  longer  or 
shorter,  depending  upon  whether  they  are  reflections  of  the 
upper  or  lower  surface,  respectively.  We  then  connect  the 
source  images  to  the  receiver  by  straight  lines.  If  a  given  line 
passes  through  all  the  surfaces,  then  the  ray  path  it  represents 
is  allowed.  If  the  line  misses  even  one  surface,  the  ray  is  not 
allowed.  In  the  figure  the  ray  AF  is  allowed.  The  other  ray 
A'F  is  not  allowed  (it  misses  at  D').  Since  the  reflection 
coefficient  2  or  .Mx  o  is  obviously  zero  at  any  surface 
which  the  ray  misses,  then  any  term  in  Eq.  (10)  in  which  that 
coefficient  appears,  e.g.,  Eq.  (Al),  will  be  zero.  Its  contribu¬ 
tion  to  the  total  pressure  field  at  the  receiver  will  also  be 
zero,  i.e.,  the  ray  does  not  exist. 
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